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Photon number states generated from a continuous-wave light source 
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Conditional preparation of photon number states from a continuous-wave nondegenerate optical 
parametric oscillator is investigated. We derive the phase space Wigner function for the output state 
conditioned on photo detection events that are not necessarily simultaneous, and we maximize its 
overlap with the desired photon number state by choosing the optimal temporal output state mode 
function. We present a detailed numerical analysis for the case of two-photon state generation from 
a parametric oscillator driven with an arbitrary intensity below threshold, and in the low intensity 
limit, we present a formalism that yields the optimal output state mode function and fidelity for 
higher photon number states. 
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I. INTRODUCTION 

The study of nonclassical states of light has provided 
a deeper understanding of quantum fluctuations and the 
role of measurements in quantum theory and it has led 
to applications in precision metrology and quantum com- 
munication. The photon number states, or Fock states, 
play a special role, because they have vanishing intensity 
fluctuations, and their interaction, e.g., with a single two- 
level atom in an optical cavity, is particularly regular. By 
injection of atoms with properly selected excitation and 
passage times through a micro maser, it is possible to 
build up a wide range of states, including number states 
of the cavity field [l|, [1J , and recent progress in extend- 
ing the photon lifetimes in microwave cavities [|[ bring 
promise for further experimental progress in this direc- 
tion. 

By control of single photo emitters such as a single 
molecule a color center j3i], or a quantum dot @, it 
is also possible to generate traveling light pulses in the 
optical frequency range that contain only a single pho- 
ton. For a review on single photon emitters see |7|. In 
these schemes, however, the production of higher num- 
ber states is not straightforward, and in the present pa- 
per we shall address an alternative conditional approach 
where the signal beam from a nondegenerate optical 
parametric oscillator (OPO) is projected into the desired 
state by a quantum measurement performed on the idler 
beam. This method has recently been used to generate 
single- and two-photon states from an OPO driven by 
a pulsed pump field Q. Conditional generation of non- 
classical states was proposed by Dakna et al @, see also 
[lOl . [ll[ , and generation of single-photon and Schrodinger 
cat states has been demonstrated in experiments where 
measurements performed on a small fraction of the light 
beam from a degenerate OPO caused the projection of 
the remaining beam 

E1H0- Single-photon and cat 
state production from continuous-wave OPOs has been 
studied theoretically in [15l.ll6l.ll7l]. and in the present pa- 
per we generalize the analysis to continuous- wave genera- 
tion of higher photon number states. In contrast to single 



photon and Schrodinger cat state production, generation 
of states with two or more photons involves multiple con- 
ditioning photo detection events, and in the continuous- 
wave case, we are particularly interested in determining 
how the generated state is affected by the temporal sep- 
aration of the conditioning detections. We analyze this 
feature generally for two-photon states and we present 
an analytical treatment, restricted to the low intensity 
limit, for n-photon state generation. 

Figure [T] exemplifies the experimental setup used to 
generate Fock states. A nondegenerate OPO produces 
pairs of distinguishable photons. The two kinds of pho- 
tons are separated to produce two correlated twin beams. 
One beam (denoted the trigger) is observed with an 
avalanche photodiode (APD) detector. Detection of n 
close detector clicks in the trigger arm projects the state 
in the signal beam into an n-photon state, and the gener- 
ation is verified by homodyne detection. Since the field is 
a continuous-wave field the temporal mode occupied by 
the produced state needs to be specified, and the largest 
n-photon fidelities are obtained by optimizing the choice 
of signal state mode function. We shall investigate to 
which extent the n photons in the signal beam may oc- 
cupy a single mode despite the click events happening 
not exactly simultaneously. 

In order not to miss close clicks due to the finite dead 
time of real detectors it might be advantageous to split 
the trigger beam and send it onto more than one APD 
detector, as was done in the two-photon experiment with 
pulsed fields [8|. However, we show that our theoretical 
expression for the conditional output state, and thus also 
the fidelity and the optimal signal state mode function, 
is independent of the number of (ideal) APD detectors 
used. 

In Sec. II we start out with a two-mode treatment of 
the two-photon state generation process. In Sec. Ill we 
generalize to the multi-mode case valid for continuous- 
wave fields. We determine the Wigner function for the 
output state conditioned on two trigger detector click 
events, and we calculate the two-photon state fidelity as 
a function of the signal state mode function. In Sec. IV 
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FIG. 1: Experimental setup for conditional preparation of 
Fock states from a type II (i.e. polarization) nondegenerate 
OPO. PBS: polarizing beam splitter, APD: avalanche photo- 
diode, BS: beam splitter, and LO: local oscillator. 



we optimize the signal state mode function over all real 
functions to obtain maximal two-photon state fidelity. 
Finally, in Sec. V, we consider n-photon state generation 
in the low intensity limit. We describe the state produced 
in the signal beam in terms of photons occupying specific 
temporal modes, and we determine the optimal output 
state mode functions. Sec. VI concludes the paper. 



II. OUTPUT STATE CONDITIONED ON TWO 
TRIGGER DETECTOR CLICK EVENTS - 
TWO-MODE TREATMENT 

In this section we describe the two-photon state gen- 
eration in the context of a simple two-mode theory to 
introduce some of the basic ideas. This treatment is ap- 
proximately valid when a pulsed pump field is used. The 
initial state generated by the nondegenerate OPO is a 
two-mode squeezed vacuum state [l8[ 
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where r is the squeezing parameter and the first (second) 
quantum number inside the ket on the right hand side is 
the number of photons in the trigger (signal) mode. 

We assume that a trigger detector click results in the 
transformation p — ► a t pa\ /Tr (fit pd\) of the density oper- 
ator p, where fit is the trigger mode annihilation operator. 
We apply the click transformation twice to the state fl}, 
and since we do not subject the trigger mode to further 
measurements, we trace over the trigger mode afterward 
and renormalize to obtain the conditional single-mode 
output state 
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The vacuum and the single-photon state contributions 
are eliminated by the conditioning procedure, and the 
generated state is a superposition of a two-photon state 
and higher photon number states. The two-photon state 
fidelity is easily obtained from ([2]) as 



(2|p c z lc fc|2) 
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The fidelity approaches unity in the limit where the 
squeezing parameter is small, because a small r corre- 
sponds to a weak pump field, and hence the probability 
to produce more than two photon pairs within a single 
pulse is small. 

For the multi-mode case it turns out to be convenient 
to describe the initial unconditional state and the condi- 
tional state in terms of Wigner functions, and we hence 
introduce this alternative approach now. Since the OPO 
is a Gaussian light source, the two-mode Wigner function 
for the initial unconditional state is a Gaussian 
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y = (xt,pt,x s ,p s ) T is a column vector of quadrature 
variables for the trigger and the signal mode, and V 
is the covariance matrix. In terms of the operators 



y = {x t ,pt,x t 
and p t/s = - 



,p s ) , defined as 



't/s 



(a t / s 



t/S 



)/V2, 



[ a t/s 



t/s 



)/\/2, where fi s is the annihi- 



lation operator for the signal mode, the elements of the 
covariance matrix are Vij — (y%yj) + (yjDi), and from (JTJ) 
we obtain 

V u = V 22 - V 33 = V M = cosh(2r), (5) 
V 13 = V 31 = -V 2i = -V^ = sinh(2r), (6) 

while the other matrix elements are zero. 

The transformation from |T]) to @ is translated into 
(see [II [l7|) 
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where N c u c k is a normalization constant. The two- 
photon state fidelity of the generated state is given by 

F 2 = 2-7T JJ W c i lc k{x s ,p s )W n=2 (x s ,p s )dx s dp s , (8) 

where W n=2 is the Wigner function for a two-photon 
state. Equation ([5]) once again leads to the result ([3]). 
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III. OUTPUT STATE CONDITIONED ON TWO 
TRIGGER DETECTOR CLICK EVENTS - 
MULTI-MODE TREATMENT 



In the continuous- wave case the field annihilation oper- 
ators are time-dependent and satisfy the commutator re- 
lation [a(t), a'(t')] = S(t — t'). In the following we denote 
the trigger beam annihilation operator a + (t) and the sig- 
nal beam annihilation operator a_ (i) to distinguish them 
from the single mode operators in the last section. In 
principle, there are now infinitely many modes, but since 
we can trace out all unobserved modes, we only need 
to consider the two trigger modes, in which the condi- 
tioning trigger detector clicks occur, and the signal mode 
occupied by the generated state, which is a great simplifi- 
cation. It is necessary to include two trigger modes since 
the conditioning clicks may happen at different times. 
Initially we assume that the trigger modes are distinct. 

The temporal shapes of the relevant modes are given 
by the mode functions /»(*), i = 1,2,3, where 1 and 2 
are trigger modes while 3 is the signal mode. The single- 
mode operators (corresponding to a t and a s ) are then 
given by 



&i = J Mt')a+{t')dt' , i= 1,2, (9) 

as = / h(f)&-(t')df. (io) 



OPO output. These are [l9( 
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(a±(t)a ± (t')) = (ai(t)a T (t r ))=0, 
where A — 2 _l c i, — 1 



(11) 



^ + e, fi — 2 — e ' 6 ^ s ^ ne nonlinear gain 
coefficient of the OPO, and 7 is the OPO output mirror 
leakage rate introduced above. Note that the dimension 
less twin beam intensity (a±(t)a±(t)) is an increasing 
function of e/7. 

In analogy to Eq. |7]) the single-mode Wigner function 
for the state conditioned on two trigger detector clicks is 



W c i lck (x 37 p 3 ) = N chck Jj\ 
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where & = D^ir, C 2 = £>i - V 55 D 2 , C 3 
W 55 V? 5 V 2 5 , C 4 = V 2 5 V 2 5 , C 5 = (V 55 )-\ and 

Di = V&[(V n - 1)^-1)+^], 

D 2 = V 55 {2V 15 V 35 (V 13 V 55 ~V 15 V 35 ) 

+^55 [V 2 5 (V 33 - 1) + V 2 5 (V U - 1)] }. 
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[oi,a\] = 1 implies that J \f t (t')\ 2 dt' = 1. We as- 
sume that the trigger modes are top hat functions of 
infinitesimal width At c and height 1 / y/ At c centered at 
the ith detection time t C i. This is valid if the duration 
of a detection is much smaller than the inverse of the 
leakage rate 7 of the OPO output mirror. The signal 
mode function is used to specify the output state, and 
hence it may be chosen arbitrarily. In Sec. IV we use 
this freedom to maximize the two-photon state fidelity. 
We assume throughout the two-photon state analysis 
that the signal mode function is real. Imperfect detec- 
tion may be taken into account by replacing a±(t) with 
y/T)t/7a±(t') + yjl — T] t / s a± A , ac (t') in equations ([9]) and 
(fT0|) . where r} t / s is the trigger/signal detector efficiency 
and d± >vac are field operators acting on vacuum. 

With the single mode operators established we proceed 
as in the second part of the last section, but since three 
modes are now included, the covariance matrix V is 6 x 
6, y = (xi,pi,x 2 ,P2,X3,P3) T , and tt 2 is replaced by n 3 
in Eq. (J4j) - To calculate the covariance matrix elements 
in terms of /3(A), rj t , rj s , and OPO parameters we need 
the two-time correlation functions for the nondegenerate 



W c iick is independent of rj t . 

If the trigger beam is divided into m beams, the field 
operator representing the field in the jth beam may be 
written as Cjoci + (t) + YliLi c jiO>i,vac(t), where Cji are 
coefficients determined from the precise arrangement of 
beam splitters, and di tVac (t) are field operators represent- 
ing vacuum states. If a click is observed in the jth and 
the fcth trigger beam in the temporal modes fi(t) and 
f-z{t), respectively, and we trace over all modes except 
the output mode (denoted by Tr ), the density operator 
Ptot is transformed into 



Ptot -* Tr' 



m — 1 



J fi(t) \Cj a + (t) + E c ji®>i,va,c{t)j dt 
j f2(t') ^C fe0 a+(t') + E C kiO-i,vac(t')^ dt' 

Ptot J W) U al(t")+ E4AWOJ dt" 



h(t"') c* al(t'") + E c*4 wc (i'") df 



(13) 
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where a normalization factor is omitted. The density 
operator p tot is the direct product of the density oper- 
ator for the OPO output and the density operators for 
the vacuum states coupled into the system via the beam 
splitters. Since the annihilation operator acting on a vac- 
uum state is zero, (fT3|) simplifies to 



Ptot 



P123 



Tr 
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h{t)c ]0 a+{t)dt / f 2 {t')c k0 a + (t')dt' 



f 2 (t")c* k0 a{(t")dt" / h{t'")c* Q al{t'")dt 



(14) 

where the trace is now over the two trigger modes, and 
pi23 is the density operator for the two trigger modes 
and the output mode. The factor CjQCkocl Q c* Q is irrele- 
vant because the transformed density operator has to be 
normalized. The conditional output state is thus inde- 
pendent of the trigger detector configuration, and it is 
justified to use the simple setup in figure Q] in a theoreti- 
cal treatment. One can verify that the Wigner function 
of the conditional state continuously approaches the out- 
come of a two-photon detection in a single trigger mode, 
when the click separation approaches zero. The argu- 
ments are immediately generalized to the case of n con- 
ditioning clicks. 

Finally, the two-photon state fidelity of the produced 
state is obtained from the conditional Wigner function 
Wdick as in Eq. © 
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The fidelity depends on the choice of signal mode function 
,fz{t). In the next section we first determine the optimal 
mode function f op (t), which leads to the largest fidelity, 
and then we present explicit results for the predictions of 
equation (fT5|) . 



IV. OPTIMAL SIGNAL MODE FUNCTION 
AND TWO-PHOTON STATE FIDELITY 



In the low intensity limit (|15[) reduces to 
limF 2 (/ 3 (i)) 



V 2 V 2 



(16) 



2[(Fii-l)(l/ 3 3-l) + ^ 1 2 3 ]' 
and for small e/7 

V 15 = J fsWe-W-^dt', (17) 

V35 = 2e^/At cVtVs J fo^e-W-^dt', (18) 
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FIG. 2: Optimized mode functions for \t C 2 — t c i\y — 4 and 
T)s = 1. The three curves correspond to e/7 = (dashed 
line), e/7 = 0.08 (solid line), and e/7 = 0.2 (dotted line). 



while V11, V33, and V13 are always independent of fs(t). 
The optimal signal mode function at very low intensity 
is thus easily obtained by optimization of V15V35 under 
the constraint J f op (t) 2 dt — 1. This leads to 



lim/ op (t) =N 



7 

cxp ( - -\t - t c i 



exp 



t c 2 



(19) 



where N is a normalization constant. 

For larger intensities the fidelity can be optimized nu- 
merically by varying the shape of the mode function un- 
til no further increase in fidelity is obtained. Optimized 
mode functions for \t C 2 — t c i\"f = 4 and different intensi- 
ties are shown in figure [H The peaks of the mode func- 
tions become sharper with increasing intensity, and a dip 
appears on each side of the peaks. This behavior is quali- 
tatively the same as what was found for the single peaked 
mode function of single-photon state generation in [17j . 

Figure [3] shows the optimized fidelity (solid lines) as a 
function of e/7 for t c \ = t C 2 and the fidelity calculated 
from the optimal mode function at zero intensity (dashed 
curves). The solid and dashed curves are almost identical 
for small e/7, so the optimal mode function at zero in- 
tensity is close to optimal in the region where the fidelity 
is high and hence provides an analytical approximation 
to the optimal choice of signal state mode function. The 
figure shows that the fidelity decreases when the inten- 
sity increases. This is as expected because a larger mean 
photon flux results in larger contributions from higher 
photon number states to the output state. 

The fidelity decreases when the temporal distance be- 
tween the trigger detector clicks increases from zero as 
is apparent from figure [H which shows the fidelity as a 
function of \t C 2 — i c i|7- When the click distance is dif- 
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FIG. 3: Fidelity for t c \ = t C 2 calculated using the optimized 
signal mode function (solid lines) and the optimal mode func- 
tion at zero intensity (|19|l (dashed lines). Perfect signal de- 
tection r/ s = 1 is assumed for the curves approaching 1 to the 
left, while g a = 0.8 for the curves approaching 0.64. 



FIG. 4: Fidelity as a function of temporal distance between 
the trigger detector clicks for e/7 = (upper solid line) and 
e/7 — 0.08 (lower solid line) calculated using the optimized 
signal mode function. The dashed line shows the fidelity for 
e/7 = 0.08 obtained from the mode function (|19|) . Perfect 
signal detection 77,, = 1 is assumed. 



ferent from zero, the two photons no longer belong to 
precisely the same mode. We return to this point in sec- 
tion V, where we also obtain an analytical expression for 
the two-photon state fidelity as a function of click dis- 
tance in the small intensity limit. The fidelity calculated 
from the optimal mode function at zero intensity (|19[) is 
also shown in figure 2] (the dashed curve), and it is seen 
that (|19p is a good approximation to the optimal mode 
function even if the click distance is large (provided e/7 
is not too large). 

We note that the generation method described here 
favors close click events since 



(a+(tci)a f + (tc2)a+(tc2)g+(tci)) _ 
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increases from one to two when \t c2 — i c i|7 decreases from 
infinity to zero: the trigger events are bunched in time. 



V. MODE OCCUPATION DESCRIPTION OF 
THE CONDITIONAL STATE IN THE LOW 
INTENSITY LIMIT 

In the previous sections we characterized the output 
state by the Wigner function W c u c k from which we were 
able to calculate the fidelity for an arbitrary state in 
an arbitrary mode. However, a deeper understanding 
of the nature of the conditionally produced signal beam 
state can be obtained by considering the state as built 



up of photons occupying specific temporal modes. A de- 
tailed mode description of multi-photon states and ma- 
nipulations of such states is given in Ref. [20(. In the 
present section we use this approach to investigate the 
state generated when conditioning on n trigger detector 
click events. It is assumed throughout that e/7 <C 1. 

If the t rigg er detector clicks are far apart, we know 
from Ref. [17 1 that the fidelity for a single photon in each 
of the n modes 



i\t-t ci \ 



1,2, 



,,n 



(21) 



is unity. On the other hand, for n = 2 and t c \ = t c2 we 
found in the last section that the fidelity is unity for two 
photons in the mode g\ (t). In both limits the state gener- 
ated in the signal beam conditioned on two clicks is thus 
on the form \tp 2 ) = N^ 2 JJ dtdt' 'g 1 (t)g 2 (t')a'_(t)a'_(t')\0) , 
where N^ 2 is a normalization constant. One is thereby 
led to consider whether this result is also valid for inter- 
mediate separation of the trigger detector clicks. In the 
appendix we show that the state generated in the sig- 
nal beam when conditioning on n trigger detector click 
events is 



n . 



|0), (22) 



where N^ n is a normalization constant. 

To illustrate the meaning of (j2"2")l we consider the case 
n = 2 in some detail. Two orthogonal mode functions 
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FIG. 5: Mode functions gi(t) and f/2(i) (solid lines) and f a (t) 
and ft(t) (dotted lines) for (tc2 — t c l)j = 8. 



are constructed from g\ (t) and g 2 (t) : 
1 



fait) = 

hit) = 



a/2(1 + /i 2 ) 
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[ffl(t) + <&(*)] 



foi(t)-0>(*)] 
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(24) 



where 



9i(t)gj(t)dt 



1 + 



t C j|J exp - <c 



(25) 



The four mode functions are illustrated in figure [5] for 
7(*c2-*ei) = 8. Inserting J3SJ) and ([M]) in we obtain 
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where \x,y) a b = \%)a ® \y)b an d \ x )j means x photons in 
the mode fj (t) . Figure [HI shows the norm square of the 
coefficients in (|2^|) . Since f a (t) is identical to the optimal 
mode function for e/7 « 1 (I19p . the upper curve shows 
the maximal two-photon state fidelity as a function of 
temporal distance between the conditioning clicks (i.e. 
the same curve as in figured]). For small 7|i c2 — i c i| we 
find from equations ([2l>]) and (f2"6")) that 



(27) 



F 2 (/^)^i-(j*c2-*ciiy 



Hence we still have a good two-photon state even if the 
trigger detector clicks are not exactly simultaneous. 

The optimal mode function for a general n is obtained 
by maximizing the fidelity between (|22p and an n-photon 
state 



1 



dtf(t)al(t) 



|0>, 



(28) 




FIG. 6: The upper (lower) curve shows the probability to 
detect two (zero) photons in the mode f a (t) and zero (two) 
photons in the mode /&(t) as a function of temporal distance 
between trigger detector clicks when e/7 — * 0. 



where f{t) is the mode function to be optimized. The 
n-photon state fidelity is 

^(/(t)) = |(^n|n) / | 2 =n!|iV J 2 n / f(U)g t iU)dt 

1=1 

(29) 

It is apparent that the phase of f(t) should be time inde- 
pendent to maximize F n (f(t)), and hence we choose f(t) 
real. Variational optimization of Yli=i I f(ti)9i(ti)dti 
leads to 

n n „ 

*/(*) = e^w n j wM^, (30) 



and thus 



(31) 



where the constants Ci and £ are determined from the 
highly nonlinear set of equations 



(32) 



j=i fe=i 



and 



1 — E] CiCjlij . 

i=l j=i 



(33) 



For n = 2 equations (|3"2"|) and (j3"3")) are £ci = ci/12 + c 2 , 
£c 2 = cx + c 2 /i 2 , and 1 = c 2 + c 2 . + 2cic 2 ii 2 , and hence 



Ci = c 2 = 1/^/2(1 + Ii 2 ) and £ = 1 + J12 in agreement 
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gives rise to the largest overlap was found. In the low in- 
tensity limit, we showed that the state generated in the 
signal beam when conditioning on n trigger detector click 



events is \ip n ) — N^ n 



|0), where 



U7=i Sdti9i(U)aUu) 

gi(t) is a function localized around the time for the ith 
click. From this we obtain the optimal mode function 
and fidelity for n-photon state generation at low inten- 
sity. For small temporal distance between the trigger 
detector clicks and n = 2 the optimized fidelity is unity 
minus a small correction of fourth order in the temporal 
click distance. 

In the present treatment we averaged over all possible 
trigger detector outcomes outside the small time windows 
specified by the trigger mode functions, but for finite 
intensities larger fidelities are obtained if we condition 
on dark intervals between the n clicks (see 



FIG. 7: Three-photon state fidelity as a function of temporal 
distance between the first and the last trigger detector clicks 
for t C 3 — t C 2 = t C 2 — t c i (solid line) and t c i ~ t C 2 (dashed 
line). The insets show the optimal mode function f(t) for the 
former case for \t C 3 — t c i|7 = 3 (left) and \t C 3 — t c i\ r y = 15 
(right), e/7 <C 1 is assumed. 



with Eq. (|19[) . For n — 3 and t C 2 — t c \ = £ C 3 — tc2 we 
obtain 



ci 



C3 = 



c 2 = -2ci/i2 4 
and for t C 2 = t c \ 



! P l2 - 2(1 + J 13 ) + 712^2 + 4(1 + hs) 
6[2/ 1 2 2 -(l + / 13 )](l+/i3) 

l + 2 C 2[2/ 1 2 2 _(i + /l3 )], (34) 



Cl 



C3 



C2 = 



(4-i? 3 ) - v/(4-/ 1 2 3 ) 2 -16(l-A 2 3 ) 



24(1 - A 2 3 ) 



l-6c 2 

3ci/i3 



(35) 



The three-photon state fidelity for these special cases is 
shown as a function of click distance in figure [7] In the 
former case the fidelity decreases from unity to 3!/3 3 , 
when the click distance increases from zero to infinity, 
while in the latter case it decreases from unity to 3-2 2 /3 3 , 
but in both cases a broad region with a fidelity close to 
unity exists. 



VI. CONCLUSION 

In conclusion we presented a theoretical description 
of conditional higher photon number state generation 
from a continuous-wave light source. We calculated the 
Wigner function for the output state conditioned on two 
trigger detector clicks and determined its overlap with a 
two-photon state. The output state mode function that 



APPENDIX A: CONDITIONAL SIGNAL BEAM 
STATE AT LOW INTENSITY 

In this appendix we prove that the state generated in 
the signal beam when conditioning on n trigger detector 
clicks at times t c \ , • • • , ten is given by (|22[) in the limit 
e/7 — > 0. Since the normally ordered moments deter- 
mine the phase space P function [2l[, it is sufficient to 
prove that the expectation value obtained from (f2"2")) of all 
normally ordered products of signal beam field operators 
is equal to the expectation value obtained from the two- 
time correlation functions (jlip , when conditioning on the 
trigger detector clicks, i.e. 



ssu[n a t(*«-)][n a -(«j)][nM*£) 

i=l j=l k=l 

n n n 

[ n «+M ) /( [n «+(*«)] [ n «+m 



q=l 



i=l 



9=1 



n n n 

[flafe)][ii (oi[n a -^) 

i—l q—1 i—1 

1 p n 

[n a -(*i)] [riM^o] [n«-(w]io) 

j=l fc=l q=l 

In „ n n 

n^)/ / [ n^i (**)][ 11 



9 =1 



2)1 



i[n«-(^)][n«-(w]io)(n^ 

i=l 



(Al) 



q=l 



The left hand side is evaluated using Wick's theorem for 
Gaussian states with zero means [22J, which states that 

(A 1 A 2 . . . A 2k ) = (iii 2 ) (A 3 A A A 5 . . . A 2k ) 
+ (A 1 A 3 )(A 2 A 4 A 5 ...A 2k ) + ... 

+ (A 1 A 2k )(A 2 A 3 A 4 ...A 2k _ 1 ), (A2) 



8 



while (A1A2 . . . A2k+i) = 0, where Ai is either an annihi- 
lation or a creation operator and k is a positive integer. 
The left hand side may thus be expressed in terms of 
(alt (f)a± (<;')) and (a±(t)a T (t')}. For small e/7 



where the summations are over all permutations of the n 
i-indices and all permutations of the n j-indices, respec- 
tively. The right hand side of (|A1|) is evaluated using the 
relation 



(aUt)a ± (t')) = ^ (l + llt-t'Qe"*!*-*'!, (A3) 



(a±(t)a T (t')) = J*£Jle-W-* 



(A4) 



It follows that the left hand side of (|Al|) is zero if m 7^ p, 
since in that case we have to combine operators, where 
the expectation value of the product of the operators is 
zero. It also follows that the denominator on the left hand 
side of (|A1[) is proportional to e 2n . In the numerator the 
lowest order terms in e are obtained by combining as 
many dl operators as possible with operators and as 
many d_ operators as possible with d + operators. These 
terms are proportional to e max ( 2 ™> 2 ™) f or m — p. To 
obtain a nonzero left hand side in the limit e/7 — * 0, 
it is thus necessary to require that m = p < n. It is 
immediately apparent that the right hand side of (|A1[) is 
zero if m — p < n is not satisfied. 

We now consider the case m = p < n and by combining 
the operators as described above and using (|A4[) and (|A3|) 
and the definitions pTjl and (|25|) we obtain for the lowest 
order terms of the left hand side (LHS) of (|A1|) 

LHS = 1 E 9h • • • 9i m (t' m )9n • • ■ 

(n — my ^— ' z — ' 



p t Pi 



9i, 



ii s -(*< 



10) 



i=l 



|0) (A6) 



repeatedly, which leads to 



01 \\a-(u) n n «-^') n «-(wiq) - 

i=l j=l fe=l g=l 

5{t" - t n+jl ) . . . 5(t m - t n+jm )5(t lm+1 - t n+jm+1 ) 

...S(t in -t n+jn ). (A7) 



Inserting (|A7[) in (|Aip we immediately obtain the result 
(|A5]l for the right hand side of ([AT]) . 
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